A Parallel Incompressible Navier-Stokes Solver with a parallel Multigrid Elliptic Kernd

We describe numerical algorithms and parallel implementations of a time-dependent, incompressible Navier-Stokes flow solver and a multigrid elliptic solver which is also used as a computation kernel in the flow solver. ~’he flow solver is based on a second-order projection method (Bell et. al [2]) applied to a slaggered finite-difference grid. The multigrid elliptic solver uses the full V-cycle scheme (13riggs, [4]) and was designed to be a general-purpose elliptic solver working on several types of finite-difference grids and for different boundary conditions. A domain-decomposition strategy is used in the parallel implementations for both the flow solver and the mtdtigrid elliptic solver on all tine and coarse grids. The implemented solvers are numerically stable and computational] y efficient, and they scale very well to a large number of processors on Intel Paragon and Cra y T3D for problems with moderate granularity. The solver codes are portable to parallel systems that support MPI, PVM and NX for interproccssor communications.